Collective stimulated Brillouin scatter 
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We develop a statistical theory of stimulated Brillouin backscatter (BSBS) of a spatially and 
temporally partially incoherent laser beam for laser fusion relevant plasma. We find a new collective 
regime of BSBS which has a much larger threshold than the classical threshold of a coherent beam in 
long-scale-length laser fusion plasma. We identify two contributions to BSBS convective instability 
increment. The first is collective with intensity threshold independent of the laser correlation time 
and controlled by difi^raction. The second is independent of diffraction, it grows with increase of 
the correlation time and does not have an intensity threshold. The instability threshold is inside 
the typical parameter region of National Ignition Facility (NIF). We also find that the bandwidth 
of KrF-laser-based fusion systems would be large enough to allow additional suppression of BSBS. 
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Inertial confinement fusion (ICF) experiments require 
propagation of intense laser light through underdense 
plasma subject to laser-plasma instabilities which can 
be deleterious for achievement of thermonuclear target 
ignition because they can cause the loss of target sym- 
metry, energy and hot electron production Among 
laser-plasma instabilities, backward stimulated Brillouin 
backscatter (BSBS) has long been considered a serious 
danger because the damping threshold of BSBS of coher- 
ent laser beams is typically several order of magnitude 
less then the required laser intensity ~ lO^^W/cm^ for 
ICF. BSBS may result in laser energy retracing its path 
to the laser optical system, possibly damaging laser com- 
ponents [H, 0. 

Theory of laser-plasma interaction (LPI) instabilities 
is well developed for coherent laser beam [3|. However, 
ICF laser beams are not coherent because temporal and 
spatial beam smoothing techniques are currently used 
to produce laser beams with short enough correlation 
time, Tc, and lengths to suppress speckle self-focusing. 
The laser intensity forms a speckle field - a random in 
space distribution of intensity with transverse correla- 
tion length Ic — FXq and longitudinal correlation length 
(speckle length) L speckle — 7F^Ao, where F is the optic 
//# and Aq = 27r/fco is the wavelength (see e.g. 0,0]). 
There is a long history of study of amplification in ran- 
dom media (see e.g [a, 0| and references there in) . For 
small laser beam correlation time Tc, the spatial insta- 
bility increment is given by a Random Phase Approxi- 
mation (RPA). Beam smoothing for ICF typically has Te 
much larger than the for the regime of RPA applicability. 
There are few examples in which the implications of laser 
beam spatial and temporal incoherence have been ana- 
lyzed for such larger Te- One exception is forward stim- 
ulated Brillouin scattering (FSBS). Although FSBS for 
a strictly coherent laser beam is a classic linear theory, 
we have obtained 0, Q its dispersion relation for laser 
beam correlation time Te too large for RPA relevance, 



but Te small enough to suppress single laser speckle in- 
stabilities . We verified our theory of this " collective" 
FSBS regime with 3D simulations. Similar simulation re- 
sults had been previously observed [nj. This naturally 
leads one to consider the possibility of a collective regime 
for BSBS backscatter (CBSBS). We will present 2D and 
3D simulation results as evidence for such a regime, and 
find agreement with a simple theory that above CBSBS 
threshold, the spatial increment for backscatter ampli- 
tude Ki, is well approximated by the sum of two con- 
tributions. The first is RPA-like cx Te without intensity 
threshold (we neglect light wave damping). The second 
has a threshold as a function of laser intensity. For Na- 
tional Ignition Facility NIF parameters the threshold is 
comparable with NIF intensities. That second contribu- 
tion is collective-like because it neglects speckle contribu- 
tions and is only weakly dependent on Te. CSBSB thresh- 
old is applicable for strong and weak acoustic damping 
coefficient Via- The theory also provides a good quan- 
titative prediction of the instability increment for small 
i^ia ^ 0.01 which is relevant for gold plasma near the wall 
of hohlraum in NIF experiments [l|. 

Assume that laser beam propagates in plasma with 
frequency ujq along z. The electric field £ is given by 
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where E{v,z,t) is the envelope of laser beam and 
B{r,z,t) is the envelope of backscattered wave, r — 
{x,y), and c.c. means complex conjugated terms. Fre- 
quency shift Aw = — 2fcoCs is determined by coupling of 
E and B through ion-acoustic wave with phase speed Cg 
and wavevector 2fco with plasma density fluctuation Sue 
given by ^ = 1(76^*'"'^+*^'^* -|- c.c, where a{Y,z,t) is 
the slow envelope and Ue is the average electron density, 
assumed small compared to critical density, Uc- The cou- 
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pling of E and B to plasma density fluctuations gives 

1 



2ko 



B=h'lla*E, (3) 
4 rir 



V ~ {dx,dy), and cr is described by the acoustic wave 
equation coupled to the pondermotive force oc £^ which 
results in the envelope equation 



K 
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= -2koE*B. (4) 



The response of the slowly varying part of Sue to the 
slowly varying part of the ponderomotive force, propor- 
tional to + responsible for self-focusing, is ne- 
glected. Via — VL/2koCs is the scaled acoustic Landau 
damping coefficient. E and B are in thermal units (see 
e.g. ^]). 

Assume that laser beam was made partially incoher- 
ent through induced spacial incoherence beam smooth- 
ing [isl which defines stochastic boundary conditions at 
z = for the spacial Fourier transform (over r) compo- 
nents E(k), of laser beam amplitude Q: 



E{k,z = 0,t) = |£;k|exp[i0k(i)], 
(expilMt) - 0k' (i')]) - '^kk' exp(-|t - t'\/T,), 
|i?k| = const, k < km] -Ek = 0, fc > km. 



(5) 



chosen as the idealized "top hat" model of NIF optics 
[l^ . Here km — fco/(2i^) and the average intensity, (/) = 
= / determines the constant. 
In linear approximation, assuming \B\ <^ \E\ so that 
only the laser beam is BSBS unstable, we can neglect 
right hand side (r.h.s.) of Eq. ([2]). The resulting linear 
equation with top hat boundary condition ([5]) has the 
exact solution as decomposition of E into Fourier series, 
E{r, z, t) = J2j with S^. oc exp [i((^k, [t - z/c) + kj ■ 

r-k2z/2fco)]. 

Figures [T] show the increment K^of the spatial growth of 
backscattered light intensity (|-Bp) oc e"^**'^ as a func- 
tion of the rescaled correlation time Tc = TckoCs/4:F'^ 
(note that definition is different by a factor 1/2F from 
the definition used for FSBS [1, obtained from the 
numerical solution of the linearized equations ©-(HI) us- 
ing operator splitting method along the characteristics 
of E and B. Here and below we use dimensionless units 
with ko/km as the unit in z direction, k^/kmCs is the 
time unit and /i = iviak^/k^, Also (...) means aver- 
aging over the statistics of laser beam fluctuations ^ 
and / is the scaled dimensionless laser intensity defined 
as / = ^^^J. Figure [1^ corresponds to the 3 -f IZ? 
simulations (three spatial coordinates and t) with the 
boundary and initial conditions ^ in the limit c — )> oo 
(i.e., setting terms in ©-(ID) to be zero. Figure [TJd 
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FIG. 1: Spatial increment Hi of CBSBS obtained from numer- 
ical simulations compared with the sum of increments kb + k,<7 
(obtained by solving ((8]) and (|10ll 'l. Parameters of simulations 
are Via ~ 0.01, F — 8. (a) 3 -I- ID simulations with Cs/c — 0, 
1 — 2 (circles) and 1 = 1 (squares). Solid and dashed lines 
show K_B for 1 = 2 and / = 1, respectively. If Ko- < then 
hb + i^a is replaced hy kb- (b) 2 + ID simulations with the 
modified boundary conditions, Cs/c = 1/500, 7 = 3 (circles) 
and 1 = 1 (squares). Error bars are also shown. Solid and 
dashed lines show hb + i^a for 1 = 3 and 1 = 1, respectively. 



shows the result of 2 + ID simulations (only 1 transverse 
spatial variable is taken into account) for the modified 
boundary condition compare with the last line in ([5]) as 
l^-kl ~ k^^^ const, k < km', E\^ = 0, k > km which is 
chosen to mimic the extra factor k in the integral over 
transverse direction of the full 3 + ID problem. In that 
case c/cg ~ 500. E.g. for T), = 0.1 we typically use 256 
transverse Fourier modes and a discrete steps Az = 0.15 
in dimensionless units with the total length of the sys- 
tem Lz = 100 and a time step At = Az/c. For each 
simulation we typically have to wait ~ 10^ time steps to 
achieve a statistical steady state and then average over 
next '^10^ time steps to find k^. 

We now relate Ki to the instability increments for (B) 
and (cr*) (we designate them kb and k^, respectively). 
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In general, growth rates of mean amplitudes only give a 
lower bound to k^. First we look for Ko-. Eq. (|3]) is linear 
in B and E which implies that B can be decomposed into 
B = X^jBkj. We approximate r.h.s. of (g]) as E*B ~ 
T,jK,Bk, so that 



= -2fco > ,K.Sk, 
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(6) 



which means that we neglect off-diagonal terms 



j j' . Since speckles of laser field arise from 



interference of different Fourier modes, j ^ j', we as- 
sociate the off-diagonal terms with speckle contribution 
to BSBS 0, [H, The neglect of off-diagonal terms 
requires that during time Tc light travels much further 
than a speckle length, L speckle ^ cTc and that Tc <C tsat, 
where tgat is the characteristic time scale at^ which BSBS 
convective gain saturates at each speckle (13] , 

Eqs. Q and ([6]) result in the closed expression 
R-fia*) = ~(kl/2){ne/ne){E* Rbb(t* E) which has the 
same form as the Bourret approximation We look 
for the solution of that expression in exponential form 
Bj,a* cx Q^i'^^+^-'^-'^t) ^ then the exponential time depen- 
dence in ([5]) allows to carry integrations in that expres- 
sion explicitly to arrive at the following relation in di- 
mensionless units 



iuj + fi + IK — {i/A)k'^ 
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where l/fc„j is the transverse unit of length and vectors 
span the entire top hat ([5]), i.e. / = J2j l-^kj 
In the continuous limit N ^ oo, sum in ([7]) is replaced 
by integral which gives for the most unstable mode k = 0: 

- iuj + 11 + in + i—I\n. — - — ^ = 0. (8) 

c c Tc 

The relation ([8|) supports the convective instability with 
the increment Ko- = Im{K) > only for / > Iconvthresh, 



where /, 



convthres 



is the convective CBSBS threshold 



given by 
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^convthresh 



(9) 



In the limit c/cg — ^ oo, the increment Kg- is independent 
of Tc which suggests that we refer to it as the collective 
instability branch. For finite but small Cg/c ^ 1 and / > 
Iconvthresh there is sharp transition of function 
of Tc from for Tc = to Tc-independent value of Ha- 
That value can be obtained analytically from (|8]) for / 
just above the threshold as follows: Ki = ii{tt/A){I — 

Iconvthresh^ / {f^I !)• 



The increment is obtained in a similar way by sta- 
tistical averaging of equation ([3]) for (B) with a* from 
equation (U) which gives 



lUJ ]- IK + t — l 

C 4: K — UJ 



(10) 



Here we neglected the contribution to kb = Im{K) from 
diffraction which gives negligible correction. Equation 
(ITOl) does not have a convective threshold (provided we 
neglect here light wave damping) while k b has near-linear 
dependence on Tc : — IJ-ITc/4: for Tc < l//i which is 
typical for RPA results. It suggests that we refer as 
the RPA-like branch of instability. 

We choose w = 0.5 in ([S]) and w = in PH]) to maximize 
Ka- and K B , respectively. Equation ([8|) also predicts abso- 
lute instability for / > fi + S^i^^ +0{fj,-^) + 0{f-^Cs/c), 
which is slightly above the coherent absolute threshold 
I — fj. but here we emphasize the convective regime. Fig- 
ures[T^ and b show that the analytical expression kb + Ko- 
is a reasonable good approximation for numerical value of 
Ki above the convective threshold © for Tc ^ 0.1 which 
is the main result of this Letter. Below this threshold the 
analytical and numerical results are in qualitative agree- 
ment at best but in that case we replace kb + Hg- by 
because k^- < in that case. 

The qualitative explanation why kb + is a sur- 
prisingly good approximation to Ki is based on the 
following argument. First imagine that B propagates 
linearly and not coupled to the fluctuations of a* , so its 
source is a*E {(t*)E in r.h.s of If {a*) cx e""^ 
grows slowly with z (i.e. if (a*) changes a little over 
the speckle length L speckle and time Tc), then so will 
at the rate 2ko- • But if the total linear response 
Rbb {I^bb is the renormalization of bare response Rbb 
due to the coupHng in r.h.s of ^) is unstable then 
its growth rate gets added to k^ in the determination 
of since in all theories which allow factorization 

of 4-point function into product of 2-point func- 
tions, (S(l)i?*(2)) = i?*?*j(l,l')^(l',2')i?*?*j*(2',2). 
Here S'(l,2) = {a* {l)a{2)){E{l)E* {2)) ~ 
{a*{l)){a{2)){E{l)E*{2)) and "1","2" etc. mean 
of all spatial and temporal arguments. 

The applicability conditions of the Bourret approxima- 
tion used in derivation of and ([TU]) in the dimension- 
less units are 



(11) 



and AujB ^ {c/cs)\kb\ as well as Awo- Here 
7o is the temporal growth rate of the spatially homoge- 
neous solution which is given by 7g = (1/4) (c/cs)/Lt/. Also 
Aujg — 1/Tc is the bandwidth for a and Aa;^ is the effec- 
tive bandwidth for B. Aujb is dominated by the diffrac- 
tion in ([3]) which gives in the dimensionless units Aw b = 
c/cs. Then (ITT|) reduces to Tc < 4/(^7) and \kb\ < 1. 
Together with the condition Tc 3> L speckle /c used in 
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the derivation of ([8]) and assuming that / — Iconvthreshi 
it gives a double inequality (77r/2)(cs/c) -C Tc ^ 7r//x 
which can be well satisfied for /i ~ 5, i.e. for v^a — 0-01 
as in gold NIF plasma but not for /i ~ 50 as in low ioniza- 
tion number Z NIF plasma. Also \hb\ < 1 implies that 
I > Iconvthresh bccause Otherwise, below that threshold, 
Kb ^ —fJ- which would contradict \kb\ < 1- All these 
conditions are satisfied for Tc ^ 1/4 for the parameters 
of Figure [T] with / = 2 or / = 3 (solid lines in Figure [1]) 
but not for 1 = 1 (dashed lines in Figure [ij. Addition- 
ally, an estimate for Tc <C tsat from the linear part of the 
theory of Ref. [13] results in the condition Tc <S 8J//Z 
which is much less restrictive than the previous condi- 
tion. These estimates are consistent with the observed 
agreement between Ki = k^- + kb smd Hi from simula- 
tions (filled circles in Figure [T|) for / above the threshold 
dH]). We conclude from Figure [T] that the apphcability 
condition for the Bourret approximation is close to the 
domain of Tc values for which Hi = k.„ + kb- 

For typical NIF parameters [H, |9i|, F = 8, ndnc = 
0.1, Aq = 351nm and Cs = 6 x 10^ cm s~^ and the elec- 
tron plasma temperature Tc ~ 5keV, we obtain from Q 
that Iconvthresh — 2.2 X lO^^W/cm^ for gold plasma with 
with i^ia — 0.01, in the range of NIF single polarization 
intensities. So we conclude that for gold NIF plasma 
/ ~ Iconvthresh whilc for low Z plasma with 0.1 
/ is well below Iconvthresh- Fig- H] shows Hi in the limit 
Cg/c = 0, Tc — !■ from simulations, analytical result Kg- 
{kb = in that limit) and the instability increment of 
the coherent laser beam Kcoherent = — (^^ — ^/)^/^/2 
(see e.g. Q). It is seen that the coherent increment sig- 
nificantly overestimates numerical increment especially 
around Iconvthresh- If we include the effect of finite Cs/c = 
1/500 and finite Tc as in Fig. then has a significant 
dependence on Tc- Current NIF SA beam smoothing de- 
sign has Tc ~ 4ps which implies T ~ 0.15. In that case 
Fig. [TJd shows that there is a significant (about 5 fold) 
change in Ki between 1 = 1 and / = 3. Similar estimate 
for KrF lasers (Ao = 248nm, F = 8, Tc = 0.7pm) gives 
Tc = 0.04 [ps] which results in a significant (40%) reduc- 
tion of Ki for J = 3 compare with above NIF estimate. 

For practical application the threshold of BSBS is often 
understood as the total gain required to amplify initial 
thermal fluctuations up to \B\'^ ~ With such defi- 

nition of threshold our results indicate that the coherent 
BSBS increment significantly overestimates Ki for prac- 
tical values of Tc as can be seen from a comparison of 
Figures [TJd andjH 

Tc in NIF can be further reduced by self-induced tem- 
poral incoherence with collective FSBS decreasing the 
correlation length with beam propagation [8, 9]. For low 
Z plasma threshold for the collective FSBS is close to 
dHI [8] . As Z increases (which can be achieved by adding 
high Z dopant), that threshold decreases below © and 
might result in an increase of the BSBS threshold. 
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FIG. 2: Ki vs. / for /i = 5.12 obtained from simulations 
(squares connected by dashed line, Cs/c = and limit Tc — > 
taken by extrapolation from Tc <C 1), analytical result Kc 
(solid curve) and coherent laser beam increment ^coherent 
(dotted curve). Upper grid corresponds to laser intensity 
in dimensional units for NIF parameters and gold plasma 
Te ~ 5keV, F = 8, n^jn^ = 0.1, Uia = 0.01, Aq = 351nm. 

In conclusion, we identified a collective threshold for 
BSBS instability of partially incoherent laser beam for 
ICF relevant plasma. Above that threshold the BSBS 
increment Ki is well approximated by the sum of the col- 
lective like increment k^ and RPA-like increment k b ■ We 
found that Ki is significantly below the BSBS increment 
Kcoherent of the Coherent laser beam. 
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